Peptide pores in lipid bilayers : voltage facilitation pleads for a revised model 
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We address the problem of antimicrobial peptides that create pores in lipid bilayers, focusing 
on voltage-temperature dependence of pore opening. Two novel experiments (voltage-clamp with 
alamethicin as an emblematic representative of these peptides and neutron reflectivity of lipid- 
monolayer at solid/water interface under electric field) serve to revise the only current theoretical 
model [2] . We introduce a general contribution of the electric field as being responsible for an 
unbalanced tension of the two bilayer leaflets and we claim that the main entropy cost of one pore 
opening is due to the corresponding "excluded-area" for lipid translation. 

PACS numbers: 05.70.Np, 82.45. Mp, 87.16. D-, 87.16.Vy 



Interaction of living cell membrane with adsorbed 
molecules and the way their uptake occurs are at the 
heart of many biological issues. Among these molecules, 
antimicrobial peptides [T] attract special attention as be- 
ing the keystone of the innate immune system of multicel- 
lular organisms. Antimicrobial peptides basically cause 
lipid-bilayer permeation by producing pores. Their uni- 
versal presence in animal and plant kingdoms, their non- 
specific broad-spectrum and their elementary structure 
let us expect their action also obeys to a widespread 
and universal physical mechanism that probably puts the 
membrane behavior in central position. 

The most accepted physical model for peptide pores 
opening is based on a tension-driven mechanism that 
can be summarized as follows [5] . Prior to pore forma- 
tion, these amphipathic peptides adsorb parallel onto the 
membrane and are supposed to increase an " internal ten- 
sion" up to a given adsorption level, beyond which they 
relax this tension by penetrating into the membrane, then 
stabilize the edge of pores that spontaneously appear in 
bilayers [3] . Although this model is the best attempt to 
formalize a widespread outlook, it is still unsatisfactory 
as it ignores two points : 1) the role of temperature and 
entropy ; 2) the role of the transverse electrical field of 
the order of 25 x 10 6 V/m experienced by living cell mem- 
branes that is known to be strong enough in some cases [1] 
to induce peptide-poration. 

In this letter, we report two novel experiments : 
1) voltage-clamp study focusing on the temperature- 
voltage dependence of pore opening with alamethicin 
as an emblematic representative of antimicrobial pep- 
tides [5] ; 2) neutron reflectivity of lipid-monolayer ad- 
sorbed at solid/water interface under electrical field. Our 
results demonstrate that the membrane behavior is cen- 
tral in this problem and the key role of its entropy. We 
support the idea that volt age- induced and zero-voltage 
peptide-poration obey the same physics : the former 
guiding us to grasp the latter. In this way, we propose a 
simple model that takes over some fundamentals of the 



tension-driven mechanism but solves some difficulties. In 
particular, we clarify the above "internal tension" as be- 
ing due to a bending energy rather than a proper mem- 
brane tension ; we introduce a general contribution of 
the electric field as being responsible for an unbalanced 
tension of the two bilaler leaflets ; and finally we show 
that for a held membrane the main entropy cost of one 
pore opening comes from the corresponding exclude-area 
for the translational entropy of lipid molecules. The im- 
pact of these two last statements could be quite general, 
not only for pore opening but also for the energetics of 
incorporation of membrane proteins in lipid bilayers. 

Voltage-clamp : Pore opening into lipid-bilayer was 
detected by voltage-clamp that consist in measuring the 
ionic current as a small electric potential is applied 
across a free-standing planar membrane, which sepa- 
rates two compartments containing KC1-1M aqueous so- 
lutions. Sample preparation and experiment setup are 
described in ref. [6| . Here, alamethicin-poration was stud- 
ied on DPPC-bilayer versus temperature above the melt- 
ing point of the lipid. This was checked by measuring at 
10 Hz the temperature dependence of the resistance and 
the capacitance of the bilayer. Both are constant above 
299 K : (4.3 ± 0.1) Gf2 and (93±l)pF, respectively. The 
latter leads to a bilayer-thickness z — 50 A, assuming the 
relative permittivity of the membrane equal to 3 [7] ■ 

In our study, alamethicin was added only to one 
side of the membrane (cis-side) at a molecular ratio 
lipid/peptide of the order of 100 (estimated from the 
area to volume ratio of the device, the membrane area 
and the amount of peptides). At a given temperature 
T, the current intensity / was recorded while the voltage 
U was alternatively set to positive and negative values 
of increasing modulus (polarity refers to cis-side). For 
each section of positive voltage, the current has a signif- 
icant non-zero value, but remains almost zero for nega- 
tive voltage. This reveals the formation of pores, which 
are induced by electrical field with the proper cis-trans 
direction, and are removed upon field inversion. Fluctua- 
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tions of current have been already analyzed [BJ. Here, we 
focus on the average value and compute for each voltage- 
section the average conductance g = I/U as a function 
of T and U. Figjl] shows a typical result. At this ionic 
strength the conductance of a single pore is typically 
1 nS [BJ , so the dashed-line in Figjl] might correspond 
to the peptide-poration transition. Note that its temper- 
ature dependence is the opposite of thermally activated 
processes such as electroporation occurring without pep- 
tide. In addition, the asymmetry regarding polarity can 
only be explained by the asymmetric peptide addition. 
This makes us confident that the dashed-line in Figjl] 
coincides with the peptide-poration transition. 
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FIG. 1. Conductance g of DPPC/alamethicin vs. tempera- 
ture T and electric potential U of the alamethicin compart- 
ment. The dashed-line has a slope equal to 6mV/K. 

The temperature dependence of the transition reveals 
the role of entropy, whereas the voltage variation shows 
an electric energy term. With the pore-free state as refer- 
ence, the Gibbs free energy per pore is G = H+QU — TS, 
with H the enthalpy gain for one pore opening, S and 
QU the entropy loss and electrical work, respectively. At 
equilibrium H + QU* = T*S, leading to : dU*/dT* = 
S/Q. Here we found S/Q ~ 6mV/K, or: 

TS/kT ~ 70 x Q/e (1) 

with k the Boltzmann constant and e the electron charge. 

Due to its rodlike a-helix structure, alamethicin carries 
a global permanent electric dipole /ip = 15 eA [4 , which 
could be viewed as responsible for voltage effects. This 
idea supposes a parallel (rather than anti-parallel) orien- 
tation of peptides forming a pore. If this were correct, as- 
suming 6 peptides per pore [5] gives Q — 6^p/z and EqJI] 
would lead to TS ~ 125 kT. This is nonphysical and in 
turn pleads for of an antiparallel orientation of dipoles so 
that their moments cancel. This is not so amazing since : 
1) although aligning one peptide dipole to the field would 
save ~ 1 kT, it would cost much more by confining paral- 
lel dipoles in repulsive interaction to form a pore; 2) the 
peptide dipole cannot explain the unequal sensitivity to 
polarity as it can favorably align with the field whatever 
its direction; 3) some peptides form voltage-facilitated 
pores despite them having no dipole [9J. So, the driv- 
ing force for voltage-induced peptide-poration necessarily 
originates from the membrane. 



Neutron reflectivity : Little is known about lipid bi- 
layers under electrical field. Structural effects were ob- 
served by infrared spectroscopy [10] or neutron reflectiv- 
ity [TT] on dried stacked assemblies of bilayers. For high 
electric field (> 10 s V/m) results suggest the alignment 
with the field of phospholipid zwiterionic head-groups. 
Could similar effects occur at lower field (i.e. compa- 
rable to the natural transmembrane field) for fully hy- 
drated head-groups ? Neutron reflectivity experiments 
have been reported on floating bilayers near a solid-water 
interface |12j. On such membranes, undulations that in- 
crease in amplitude with electrical field, dominate the 
reflectivity spectrum and likely hide more subtle changes. 

To overcome this issue, we performed neutron re- 
flectivity (EROS/LLB spectrometer) on single mono- 
layers adsorbed at the interface between saline heavy 
water ([KC1]=1M) and conductive silicon wafer allow- 
ing the electric field to be applied. The experimental 
setup is described in ref. [13] . Si- wafers were silanized 
following ref. [T3] using octadecylthrichlorosilane (OTS). 
DPPC-monolayer was deposited on silanized wafer in the 
"liquid-expanded" (LE) phase by a modified Langmuir- 
Shaefer technique avoiding contact of the deposited layer 
with air. Deposits were done at 25° C at the controlled 
number density of 1/86 A 2 i.e. comparable to that of 
aqueous bilayers and ensuring the fluidity of the film. 
The overall capacitance was measured as C = 1 fiF. The 
electric field in the monolayer is E = (C/Ae m )U, with 
e m = 3eo the membrane permittivity, eo that of vacuum 
and A = 12.6 cm 2 the area. One gets : E = 30 x 10 6 V/m 
for U — IV that is comparable to natural transmem- 
brane field. Reflectivity measurements were done at 
T = 23°C for wafer potential U = 0,-1, +1V. Their 
reliability was checked at least three times. We always 
observed the following behavior (Figj2| : spectra divided 
by the Fresnel reflectivity show characteristic oscillations, 
which arc fully superimposable for U — and —IV, but 
are unambiguously shifted to low values of transfer vector 
q for +1 V. This shift is irreversible at this temperature. 
Fitting in detail the reflectivity-data is out of the scope 
of this letter. Here, it is sufficient to note that the shift to 
low-g corresponds to thickening of the overall layer by : 

Az ~ (l/q* +lv - l/q* ov ) x 2tt = (3.6 ± 0.2) A (2) 

with q* the positions of first maxima. This thickening 
is the opposite of what is expected for capacitive di- 
electric compression. The asymmetry regarding polarity 
necessarily involves the only permanent dipoles in the 
system, i.e. those of head-group zwitterions of moment 
= 4eA [15] . Given the size of head-groups [IB] , Eqj2] 
is fully consistent with their orientation almost paral- 
lel to the plane at or -IV and aligned to the field for 
+1 V. Due to incompressibility, thickening goes with area 
shrinking of the monolayer. Transposed to a bilayer this 
shrinking only concerns the leaflet on the side of negative 
potential and tends to bend the membrane with the con- 
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vexity on positive-side, i.e. the side of voltage-facilitated 
peptide-poration (Figjlj). In case of a held bilayer, the 
shrinking is not real but only rises the tension of the 
corresponding leaflet and the spontaneous curvature. In 
case of a laterally free monolayer, it is real and proba- 
bly causes the transition from LE to a more condensed 
phase. This likely explains the observed irreversibility. 



Wafer polarity 




FIG. 2. Reflectivity of DPPC-monolayer at D 2 0/OTS-Si- 
wafer interface divided by the Fresnel's reflectivity of D2O/S1 
vs. transfer vector q. Dashed lines mark oscillations maxima. 

Model for peptide-poration - 1) Spontaneous pores : 
Spontaneous opening of one pore in lipid bilayer involves 
a mechanical energy [3] : 
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7 e 27iT - a m0 nr 



(3) 



with r the radius of the pore, 7 e its edge-energy and 
o~ m the bare-membrane tension. E m reaches a maximum 
E* n = 7T7e/c TO o for r* m = j e /o- m o below which the pore 
tends to close and beyond which it grows indefinitely. 

2) Tension- driven mechanism [5] : In this now classi- 
cal model, when a peptide adsorbs onto the bilayer, it 
pushes away the lipid head-groups in order that its hy- 
drophobic part meets the heart of the bilayer. Thus, a 
symmetric adsorption on both leaflets expands the area 
A of the bilayer. For small area number density x of 
peptides : AA/A = xAp/2, with Ap the area per pep- 
tide. The key point is that this is supposed to equally 
increase the membrane tension o~ m = k s (AA/A), with k s 
the stretching modulus (implicitly in ref . [2] , this amounts 
to considering the membrane as a canonical ensemble) : 



da m /dx = k s Ap/2 



(4) 



<r m would increase until E^ is small enough to allow 
thermal fluctuations to form pores larger than r* m . Let 
us denote a;* the corresponding rr-value. Assuming that 
peptides incorporated into the membrane at the pore 
edge have no effect on surface tension, it can be shown 
that further increase of x enriches the incorporated- 
peptide population only (i.e. their area number den- 
sity Xi increases whereas X X 1 - — - X IS constant), sim- 
ilarly to a phase transition. Thus, beyond x* the con- 
tribution of adsorbed-peptides to E m in Eqj3] is E ap — 
—wr 2 (x*k s Ap/2). If the line density p of peptides on the 
pore rim is constant : Xi = nip2irr, with rij the area 



number density of pores, then E ap splits into two terms : 

E ap = - (xk s Ap/2) nr 2 + (2TT 2 n lP k s Ap/2) r 3 (5) 

With Eq|3]in mind, the term in r 3 allows pores of radius 
larger than r* m to be stable. This model faces two prob- 
lems : 1) these peptides are amphipathic, so their sur- 
face excess is positive (they populate the surface rather 
than the membrane bulk) and from the Gibbs adsorp- 
tion isotherm da m /dx should be negative. Indeed, EqfJ] 
should be revised as the grand canonical ensemble is 
likely to be more adequate ; 2) Electric field applies a 
capacitive pressure on bilayer that lowers its effective ten- 
sion [17] . So, how to account for voltage- induced pores ? 

3) Spontaneaous curvature : Here, we propose to solve 
the above issues by noting that when peptides are asym- 
metrically adsorbed on one given bilayer leaflet, they 
tend to increase its "natural" area and thus the spon- 
taneous bilayer curvature, cq, with the same convexity 
orientation as cis-trans electric field does. The curva- 
ture energy stored by a flat bilayer is E e i — A\k c c\, 
with k c the bending elastic modulus. For spherical cur- 
vature Co <C Co = AA/A2z, with AA the area dif- 
ference between leaflets. For x peptides per unit area 
adsorbed on one given leaflet, we obtain : AA/A — xAp. 
Thus c = xC and E el — A\k c C?x 2 , with C = A P /2z. 
On the contrary, peptides in pores equally contribute to 
both leaflets and do not affect E e i (as for a m within the 
tension-driven mechanism). So, with one pore involv- 
ing p2irr peptides, the curvature energy is reduced to 
E el = A\k c L 2 (x - p2nr/A) 2 . The difference is : 



E n 



-k c C xp x 27rr, 



(6) 



instead of Eq(5] Note that by establishing bridges be- 
tween the two bilayer leaflets, pores increase their cou- 
pling that in turn increases the bending rigidity. This 
effect has been observed by neutron spin echo [18] . It 
can be the origin of cooperativity of pore opening. 

4) Pore edge : Although antimicrobial peptides are 
in many cases rodlike, they do not display the symme- 
try of a solid of revolution with respect to their am- 
phiphilicity. Instead, the projection in the plane nor- 
mal to their axis shows a bipolar repartition |19j . When 
peptides are parallel to the pore axis, their interfacial 
energy £p is minimized when they are at the pore edge 
with their hydrophilic zone facing the channel. Consider 
the projection of the channel in the membrane plane 
as a polygon of angle a with peptides as vertices, we 
note a* the angle of the same vertex that includes the 
hydrophilic zone of the peptide (a* < tt otherwise hy- 
drophilicity dominates and peptides likely do not ad- 
sorb). Ep is minimum for a = a* , occurring for the opti- 
mal pore radius r* — 1/ ( P (tt — a*)). Deviation from a* 
increases £p = rpzAaA-f, with rp the radius of peptide- 
rod, Aa = \\a — a*||, and A7 > the surface tension 
difference between heterophilic and homophilic contacts 
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(A7q<q* = lw-lw and A7 Q » <Q = 7™ -7^; superscripts 
h and w tag for peptide hydrophobic h or "waterphilic" w 
zones; subscripts for the facing medium i.e. either hydro- 
carbonated lipid tails h or water w). For sake of simplifi- 
cation let us assume that 7^ = 7™ = and 7^ = 7™ = 7. 
The interfacial energy per pore E- m t — p2nr£p is thus : 



E„ 



2nrpz x 7 x ||1 — {r/r*)\\ 



(7) 



Ei nt adds to the line tension in Eq[3] Doing so, it in- 
creases both and E^ , but also introduces a minimum 
for E m at r* < r* n , allowing the pore to be stable without 
needing to cross the energy barrier E^. 

5) Electrical work : The electrical field tends to orien- 
tate lipid head-groups of the "negative-side" leaflet. The 
corresponding potential electrical energy is exactly bal- 
anced by an asymmetric tension of this leaflet of energy 
^loaf = (A/Ai^niJJ cos(8)/z, with A\, the area per lipid 
and 9 the angle of head-groups dipoles with the field. 
When a pore opens, this tension is relaxed by E\ ea fKr 2 / A. 
This leads to a simple expression for the electrical work 
of lipid dipoles coming with one pore opening : 



(ttt-Vl cos(6)/A h z) U 



(8) 



Wl amounts to removing all the lipid dipoles in the area 
irr 2 . Its analogous to take off the charges accumulated 
at the surface of the dielectric, — U 2 e m irr 2 /2z, is always 
negligible compared to Wl for usual voltages. Eq|8] is 
thus the dominant electrical contribution. Note that this 
could be also relevant for peptide-free electroporation. 

6) Entropy : The translational entropy of Ni, lipids 
in an area A is kN\ J hi{eA/Ni J Ai J ). For a held mem- 
brane under tension, the opening of a pore of area irr 2 
only relaxes a bit the tension and lets the overall area 
unchanged. So, the accessible area is reduced by the 
"excluded- area" of the pore. The entropy cost is : 



TSl = kTnr 2 /A L 



As for the translational entropy of adsorbed peptides, 
it is also reduced by this excluded-area but also by the 
immobilization of peptides in pores assumed immobile. 
The resulting cost per pore is : —TSp = kTx (x/x*) irr 2 . 
As x/x* ~ 1 and x <C 1/Al, thus Sp <C 5l- It can be 
checked that it will be the same for all contributions of 
peptides to the entropy loss (rotational or conformational 
freedom etc) . This is due to lipids outnumbering peptides 
and to the extensiveness of entropy. Eq|9] is thus the 
dominant contribution to the entropy cost of one pore. 

Finally, the sum of Eq[3] and Eq(6] to [9] estimates 
the main contributions in the free energy per pore. 
In particular, this model allows us to account for the 
voltage-temperature dependence of pore opening : from 
Eq|8] and [9] one expects at the transition dU* /dT* = 
ft/O L cos(0)/z). With 9 = 77.4° [20], one obtains 



dU* /dT* = 5mV/K, in very good agreement with the 
result in Fig[l] Here, one understands that beyond an 
epiphenomenon, voltage effects amount to put the pore- 
opening transition in the correct temperature-window. 

Voltage-induced peptide pores always open with cis- 
to-trans electric field regarding peptide. In this letter, 
consistently with our experiments, we propose that this 
asymmetry, as well as the electric work, originate from 
the membrane rather than from the peptide. We also ar- 
gue that the main contribution to the entropy cost comes 
an "excluded-area" effect on lipid entropy. Thus, our re- 
sults place the membrane behavior in central position for 
the energetics of pore opening and thus could have a more 
general impact in modeling large molecule incorporation. 
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